Draft version June 28, 2011 

Preprint typeset using I^T^^X style emulatcapj v. 8/13/10 



O 
(N 

C3 

1—5 

(N 

6 



> 
oo 

o 



IMPACTS OF COLLECTIVE NEUTRINO OSCILLATIONS ON CORE-COLLAPSE SUPERNOVA EXPLOSIONS 
YuDAi Suwa\ Kei Kotake^"', Tomoya Takiwaki', Matthias Liebendorfer^, and Katsuhiko Sato'''^ 

Draft version June 28, 2011 

ABSTRACT 

By performing a series of one- and two-dimensional 2D) hydrodynamic simulations with spec- 
tral neutrino transport, we study possible impacts of collective neutrino oscillations on the dynamics 
of core-collapse supernovae. To model the spectral swapping which is one of the possible outcome 
of the collective neutrino oscillations, we parametrize the onset time when the spectral swap begins, 
the radius where the spectral swap occurs, and the threshold energy above which the spectral inter- 
change between heavy-lepton neutrinos and electron/anti-electron neutrinos takes place, respectively. 
By doing so, we systematically study how the neutrino heating enhanced by the spectral swapping 
could affect the shock evolution as well as the matter ejection. We also investigate the progenitor 
dependence using a suite of progenitor models (13, 15, 20, and 25 Mq). We find that there is a 
critical heating rate induced by the spectral swapping to trigger explosions, which significantly differs 
between the progenitors. The critical heating rate is generally smaller for 2D than ID due to the 
multidimensionality that enhances the neutrino heating efficiency. For the progenitors employed in 
this paper, the final remnant masses are estimated to range in 1.1-1.5Mq. For our 2D model of the 
15M0 progenitor, we find a set of the oscillation parameters that could account for strong supernova 
explosions 10^^ erg), simultaneously leaving behind the remnant mass close to ^ IAMq. 

Subject headings: hydrodynamics — neutrinos — radiative transfer — supernovae: general 



1. INTRODUCTION 

Although the explosion mechanism of core-collapse 
supernovae is not completely understood yet, current 
multi-dimensional (multi-D) simulations based on re- 
fined numerical models show several promising scenar- 
ios. Among the candidates are the neutrino heat- 
ing mechanism aided by convection a nd standing ac- 
cretion shock i nstability ( SASI) (e.g.. iMarek &: Jankal 
[200l iBruenn e t al. 2009; Suwa et al.l l2010D. the acous- 
tic mechanism (|Burrows et al .Il2007bt) or the magnetohy- 
drodvnamic (MHD) mechanism ( e .g..lKotake et al.ll2004 , 
200l lObergaulinger et all [200a iBurrows et all 12007a : 
Takiwaki et al.ll2009D . Probablv the best-studied one is 
the neutrino heating mechanism, whose bas ic concept 
was first proposed bv 'Colgate fc White! ()1966D . and later 
reinforced bytBethe & Wilson (1985) to take a currently 
prevailing delayed form. 

An important lesson from the multi-D simula- 
tions mentioned above is that hydro dynamic mo- 



|2010( ) can help the onset of the neutrino-driven ex- 
plosion, which otherwise fails ge nerally in spheri- 
cally symmetric (ID) s imula ti ons ([Liebendorfer et al] 
200lt iRampp fc Jankal [2002t [Thompson et all 120031: 



Sumivoshi et al.ll2005() . This is mainly because the ac- 
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cretion timescale of matter in the gain region can be 
longer than in the ID ease, which enhances the strength 
of neutrino-matter coupling there. 

In fact, the neutrino-driven explosions have been ob- 
tained in the following state-of-the-art two-dimensional 
(2D) simulations. Using the MuDBaTH code which in- 
cludes one of the best available neutrino transfer ap- 
proximations, iBuras .et_al., (20QJ) firstly reported ex- 
plosions foi_a^ji£n-rotatiiigk)w-mass (11.2Mq) progen- 
itor of IWooslev et aD (I2002D. a nd then for a 15Mq 
progenitor of IWooslev fc Weaverl (I1995D with a m oder- 
ately rapid rotation imposed ( Marek fc Ja nka"2009i) . By 
implementing a multi-group flux-li mited d iffusio n algo - 
rithm to the CHIMER A code (e.g.. IBruenn et al.ll2009| ). 
lYakunin et al.l ()2010[ ) obtained explosio ns for a non- 
rotati ng 12M0 and 25 Mf:^ progenit o r of IWooslev et al.l 
()2002D . More recently, iSuwa et al.l (I2010D pointed out 
that a stronger explosion is o btained for a rapidly 
rotati ng 13M0 progenitor of iNomoto fc Hashimotol 
(jl988l ) compared to the corresponding non-rotating 
model, in which the isotropic diffusion source approx- 
imation (IDSA) for the spectral neutrino transport 
(jLiebendorfer et aT1l2009() is implemented in the ZEUS 
code. 

However, this success opens further new questions. 
First of all, the explosion energies obtained in these sim- 
ulations are typically underpowered by one or two orders 
of magnitudes to explain the canonical supernova kinetic 
energy (~ 10^^ erg). Moreover, the softer nuclear equa- 
tion o f state (EOS), such as of the iLattimer fc Swestvl 
()1991D (LS) EOS with an incompressibility K = 180 MeV 
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at nuclear densities is employed in those simulations. On 
top of evidence that favo rs a stifFer EOS bas ed on nu- 
clear experimental data (jShlomo et al.l I2006D , the soft 
EOS may not account fo r the recently observe d massive 
neutron star of ~ 2Mq (jPemorest et al.||2010f) fsee the 
maximum mass for the LS180 EOS in " lO'Connor fc Ottl 
l20Tl[) . With a stiffer EOS, the explosion ene rgy m ay be 
even lower as inferred from iMarek k, Jankal ()2009() who 
did not obtain the neutrino-driven explosion for their 
model with K = 263 MeV. What is then missing fur- 
thermo re? We may get the an swer by going to 3D simu- 
lations (jNordhaus et alJl201Cl[ ) or by taking into account 
new ingredients, s uch as exotic physi cs in the core of the 
protoneutron star (iSagert et al.ll2009D. viscous heating by 
the magne torotation al instabilitv (jThompson et "aLll2005l : 
[Masada ct alJ 1201 or energv dissipation via Alfven 
waves (jSuzuki et alJl2008| ). 

Joining in these efforts, we explore in this study the 
possible impacts of collective neutrino oscillations on en- 
ergizing the neutrino-driven explosions. The collective 
neutrino oscillations, i.e. neutrinos of all energies that 
oscillate almost in phase, are attracting great attention, 
because they can induce dramatic obser vable effects such 
as a spectral split or swap (e.g.,|R affelt fc Smirnovll2007l : 
iDuan et"al] l2008t iDasgupta et al., .2008i . and see refer- 
ences therein) . They are predicted to e merge as a d i stinct 
feature in their energy spectra (see iDuan et "all 120101 : 
iDasguptal IMol for reviews of the rapidly growing re- 
search field and collective references therein). Among 
a number of important effects possibly created by the 
self-interaction, we choose to consider the effect of spec- 
tral splits between electron- (z^e), anti-electron neutrinos 
{ve), and heavy lepton neutrinos (vx, i.e., v^, Vj and their 
anti-p articles) above a threshold energy (e.g.. lFogli et al.l 
(|2007t l). Since v^'A have higher average energies than the 
other species in the postbounce phase, the neutrino fla- 
vor mixing would increase the effective energies of I'e and 
j^e, and hence increase the neutrino heating rates in the 
gain region. A formalism to treat the neutrino o scillation 
in the Boltzmann neutrino trans port is given in |Yamada| 
(|2Q00( ): iStrack fc Burrow s' (2005), but difficuh to imple- 
ment. To just mimic the effects in this study, we perform 
the spectral swap by hand as a first step. By changing 
the average neutrino energy, (ei/^), as well as the position 
of the neutrino spheres (i?iy^) in a parametric manner, 
we hope to constrain the parameter regions spanned by 
(e^^) and in which the additional heating given by 
the collective neutrino oscillations could have impacts on 
the explosion dynamics. Our strategy is as follows. By 
performing a number of ID simulations, we will firstly 
constrain the parameter regions to some extent. Here we 
also investigate the progenitor dependence using a suite 
of progenitor models (13, 15, 20, and 25 Af©). After 
squeezing the condition in the ID computations, we in- 
clude the flavor conversions in 2D simulations to see their 
impacts on the dynamics, and we also discuss how the 
critical condition for the collective effects in ID can be 
subject to change in 2D. 

The paper opens with descriptions of the initial mod- 
els and the numerical methods focusing how to model the 
collective neutrino oscillations (Section 2). The main re- 
sults arc shown in Section 3. We summarize our results 
and discuss their implications in Section 4. 



2. NUMERICAL METHODS 

2.1. Hydrodynamics 

The employed numerical methods are essentiahy th e 
same as those in our previous paper (jSuwa et al.ll2010[) . 
For later convenience, we briefly summarize them in the 
following. The basic evolution equations are written as, 

l + pV.v.O, (1) 

P^ = -VP-pV<i>, (2) 

de* 

— +V-[(e* + P)v] = -pvV$ + QB, (3) 
at 

f-Q.. (4) 

A $ = AirGp, (5) 

where p, v, P, v, e* , $, are density, fluid velocity, gas pres- 
sure including the radiation pressure of neutrinos, total 
energy density, gravitational potential, respectively. The 
time derivatives are Lagrangian. As for th e hydro solver, 
we employ the ZEUS-2D code (jStone fc No rman 1991 
whic h has been modifie d for core-collapse simulation s 
(e.g.. ISuwa eFall I2007blial . [2009t ITakiwaki et alllMl . 
Qe and Qn (in Equations ([3]) and (|4])) represent the 
change of energy and electron fraction (Ye) due to the 
interactions with neutrinos. To estimate these quanti- 
ties, we implement spectral neutrino transport using the 
isotropic diffusion source approximation (IDSA) scheme 
(jLiebendorfer et al.l[2009( ). The IDSA scheme splits the 
neutrino distribution into two components, both of which 
are solved using separate numerical techniques. We ap- 
ply the so-called ray-by-ray approach in which the neu- 
trino transport is solved along a given radial direction as- 
suming that the hydrodynamic medium for the direction 
is spherically symmetric. Although the current IDSA 
scheme does not yet include i^x and the inelastic neu- 
trino scattering with electrons, these simplifications save 
a significant amount of computationa l time compared to 
the ca nonical Boltzmann solvers (see iLiebendorfer et ahl 
(l2009h for m ore de tails). Following the prescription in 
IMiiller et al.l (l20Tol) . we improve the accuracy of the to- 
tal energy conservation by using a conservation form in 
equation ([3]) , instead of solving the evolution of internal 
energy as originally designed in the ZEUS code. Numer- 
ical tests are presented in Appendix A. 

The simulations are performed on a grid of 300 log- 
arithmically spaced radial zones from the center up to 
5000 km and 128 equidistant angular zones covering 
< 9 < TT for two-dimensional simulations. For the spec- 
tral transport, we use 20 logarithmically spaced energy 
bins reaching from 3 to 300 MeV. 

2.2. Spectral swapping 

As mentioned in 511 we introduce a spectral inter- 
change from heavy- lepton neutrinos (I'p, Vt and their 
antineutrinos, collectively referred as hereafter) to 
electron-type neutrinos and antineutrinos, namely i^x ^ 
Ue and Vx — > Ve- Instead of solving the transport 
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equations for v^, we employ the so-called light-bulb ap- 
proximation and focus on the o ptically thin region out- 
side the neutrinoso phere (e.g., iJanka fc Mueller! 119961 : 

lOhnishi et an[20nfih. 

According to iDuan et al.l (j2010t) . we set the threshold 
energy, eth, to be 9 MeV, above which the spectral swap 
takes place. Below the threshold, the neutrino heating is 
estimated by the spectral transport via the IDSA scheme. 
Above the threshold, the heating rate is replaced by 

/■oc 

Qe^ de, + Mr, e.), (6) 

where j and x are the neutrino emissivity and absorp- 
tivity, respectively, and fi,{r, e^) corresponds to the neu- 
trino distribution function for with being energies 
of electron neutrinos and antineutrinos. In the light-bulb 
approach, it is often approximated by the Fermi-Dirac 
distribution with a vanishing chemical potential (e.g., 
lOhnishi et al. 2006) as, 

where k, T,^^ are the Boltzmann constatn and the neu- 
trino temperature, respectively. g{r) is the geometric 

1/2 

factor, g{r) = 1 — [l — {Rv^/r)^] which is taken into 
account for the normalization, with R^^ being the radius 
of the ncutrinosplicre. The neutrino luminosity of Vx at 
the infinity is the given as 

where {e^J = de^^elJ^{e^J)/ de^^elj^{t^^)^ is 
the average energy of emitted neutrinos. The position 
where the spectral swapping sets in is fixed at 100 km 
(around the gain radius) and the onset time is varied as 
a parameter, =100, 200, and 300 ms after bounce. 

In fact, the threshold energy depends on the neu- 
trino lu minosities, spectra and oscillation parameters 
(see, e.g.. lDuan et al.ll2010l and references therein) with 
conserved net Ve flux (i.e., the lepton number conserva- 
tion). However, the conservation of lepton number is too 
complicated to satisfy in the dynamical simulation be- 
cause the neutrino spectrum and the luminosity evolve 
with time. In order to focus on the hydrodynamic fea- 
tures affected by the spectral modulation induced by the 
swapping, we simplify just a single threshold energy in 
this work. 

To summarize, the parameters that we use to mimic 
the spectral swapping are the following three items, (i) 
R^^ which is the radius of the neutrinosphere of v^, (ii) 
(cjy^ ) which is the average energy of , and (iii) ts which 
is the time when the spectral swapping sets in. 

3. RESULT 
3.1. One- dimensional models 
3.1.1. ID without spectral swapping 

In this subsection, we first outline the ID collapse dy- 
namics wit hout spectral swapping. W e take a 13 Mq 
progenitor (jNomoto fc Hashimotolll988D as a reference. 

At around 112 ms after the onset of gravitational col- 
lapse, the bounce shock forms at a radius of ~ 10 km 




Time after Bounce [ms] 



Figure 1. Time evolution of ttie neutrino luminosity (top panel) 
and average energy (bottom panel) for Ue (red-solid line) and 
(blue-dashed line). 

with an enclosed mass of - Q.7MM The central den- 
sity at this time is pc — 3.6 x 10^ g cm~'^. The shock 
propagates outwards but finally stalls at a radius of 
100 km. Due to the decreasing accretion rate through 
the stalled shock, the shock can be still pushed outward. 
However, after some time, the shock radius begins to 
shrink. The ratio of the advection timescale, Tadv, and 
the heating timescale, Theat, is an importa nt indicator 



for the criteria of neutrino driven explosion 


Buras et al.l 


l2006HMarek fc Jankall2009HSuwa et al.ll201C 


). In our ID 



simulations, Tadv/Theat is generally smaller than unity in 
the postbounce phase. This is the reason why our ID 
simulations do not yield a delayed explosion. This also 
the case for the other progenitors (15, 20, and 25 Mq) in- 
vestigated in this study. As for the accretion phase (later 
than ~ 50 ms after the bounce), the typical neutrino lu- 
minosity at r = 5000 km is 3 x 10^^ erg s~^ for both 
Ve and j/g, and the typical average energy is {e^^) ~ 9 
MeV and (epj « 12 MeV as shown in Figure[Tl Figured] 
indicates the resultant neutrino luminosity spectrum at 
100 ms after the bounce. 

3.1.2. ID with spectral swapping 

The investigated models with the spectral swapping 
are summarized in Table [TJ As already mentioned, the 
model parameters arc the neutrinosphere radius (i?j/^), 
the average energy of neutrinos {{f-vj)), and the onset 
time of the spectral swapping {ts). The model names in- 
clude these parameters; "NH13" represents the progen- 
itor model, "R.." represents R,y^ in units of km, "E.." 
represents (ci/^) in MeV, "T.." represents ts in ms, and 
the last letter "S" represents ID (spherical symmetry). 

Figure [3] presents the time evolution of the mass shells 
for models NH13R30E12T100S and NH13R30E13T100S. 
The difference between these panels is the average en- 
ergies of neutrinos, (ci/^.) = 12 MeV for the top panel 
and 13 MeV for the bottom panel. The thick solid lines 

^ Note that 0.7Mq is rather high value that is due to ap- 
proximations employed in our simulation. We omit the electron 
scattering by neutrinos and general-relat ivistic effects, which lead 
smaller inner core ma ss at bounce (see ILiebendorfer et al.l 120011 : 
[Thompson et ani2003I V In addition, mo re improved electron cap- 
ture treatment would lead even smaller ijLanganke et al]|2003l l. 
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Figure 2. The neutrino luminosity spectrum of Ve (red-solid line) 
and Pe (blue-dashed lino) without the spectral swapping at 100 ms 
after the bounce. For comparison, we show the injected luminosity 
spectrum of v^. with (e^^ ) = 15MeV, which will be swapped with 
the original spectrum of and Pe at > 9 MeV for models 
including spectral swapping. 

represent the radial position of shock waves. Regardless 
of a small difference of (e^^J, model NH13R30E13T100S 
shows a shock expansion after the manual spectral swap- 
ping is switched on (see the thick line in the bottom 
panel of Figure [3]) , while the stalled shock does not re- 
vive for model NH13R30E12T100S (top panel). This 
suggests that there is a critical condition for the suc- 
cessful explosion induced by the spectral swapping. In 
the bottom panel, the regions enclosing the mass of 
Mr ^ \.2Mq (thin black line) corresponds to the so- 
called mass cut, which could be interpreted as the final 
mass of the remnant. The fact that a clear mass cut 
emerges in model NH13R30E13T100S indicates that a 
neutron star will be left behind in t his model. Such a def - 
inite mass-cut has been observed in lKitaura et al.l (|2006f ) 
who reported a successful neutrino-driven explosion (in 
ID) for a lighter progenitor star, which is, however, diffi- 
cult to realize for more ma ssive stars in 2D ( e.g., Figure 
2 in iMarek fc Jankal (120091 ) and Figure 1 in Isiawa et al.l 

(|2om 

As a tool to measure the strength of an explosion, we 
define a diagnostic energy that refers to 



E, 



diag 



dV 



1 



(9) 



where e is internal energy, D represents the domain in 
which the integrand is positive. Figure U shows the time 
evolution of E'diag for some selected models. The diag- 
nostic energy increases with time for the green-dotted 
line, which turns to decrease for the red line, noting that 
the difference between the pair of models is A (ci, ) = 1 
MeV. The blue-dashed line (model NH13R30E15f IOCS) 
has {e^J = 15 MeV and reaches larger i^diag than the 
green line (NH13R30E13T100S; {e^J = 13 MeV). On the 
other hand, the later injection of the spectral swapping 
leads to smaller -Ediag, i-e. the brown-dot-dashed line 
{ts=200 ms) shows smaller i^diag than the blue-dashed 
line {ts =100 ms). For models that experience earlier 
spectral swapping with higher neutrino energy, the diag- 
nostic energy becomes higher in an earlier stage, as it is 
expected. 
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Figure 3. Time evolution of mass shells for NH13R30E12T100S 
(top) and NH13R30E13T100S (bottom). The Black thin line cor- 
responds to 1.2Mq and the black thick line represents the shock 
wave position, respectively. The difference between these panels 
is the average energies of neutrinos, (ei/^,) = 12 MeV for the top 
panel and 13 MeV for the bottom panel. 

Looking at Figure H] again, i?diag for the exploding 
models seems to show a saturation with time. These 
curves can be fitted by the following function. 



i?diag(t) = i?dTag(l-e 



-at+b 



(10) 



where E^^^ is a converging value of iJdiag, a and b 
are the fitting parameters. As for NH13R30E13T100S, 
^dikg = 10^° erg. This fitting formula allows us 

to estimate the final diagnostic energy especially for the 
strongly exploding models whose diagnostic energy we 
cannot estimate in principle because the shock goes be- 
yond the computational domains (r < 5000 km) before 
the saturation. 

Figure [5] shows the summary of ID models. For a given 
neutrino luminosity that is determined by R^^ and (e^^) 
(equation jS])). The gray lines correspond to the neutrino 
luminosities determined by the pairs of R,^^ and (e^^) 
which is 1 to 5 X 10^^ erg s~^ from bottom to top lines. 
Circles and crosses correspond to the exploding and non- 
exploding models, respectively. Not surprisingly, explo- 
sions are more easier to be obtained for higher neutrino 
luminosity. 

As is well known, the combination of (cj/^) and L^^ is 
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Figure 4. Diagnostic energies as functions of time. Red- 
solid, green-dotted, bluc-dashcd and brown-dot-daslicd lines 
correspond to models NH13R30E12T100S, NH13R30E13T100S, 
NH13R30E15T100S, and NH13R30E15T200S, respectively. The 
average energies of three lines other than brown-dot-dashed line 
differ from each other, i.e. (t^^) = 12 MeV (red solid), 13 MeV 
(green dotted), and 15 MeV (blue dashed), respectively. As for 
the brown-dot-dashcd line, only the onset time of spectral swap- 
ping is different from the blue-dashed line. The red-solid line shows 
only an oscillation, while the other lines show increasing diagnostic 
energy. 
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Figure 5. Summary of ID models with ts = 100 ms after the 
bounce. Circles indicate the exploding models, while crosses show 
non-exploding models. Gray solid lines correspond to the luminos- 
ity of calculated by Eq. JHJ, which are 1, 2, 3, 4, and 5 X 10^^ 
erg s~^ from bottom to top. 

an important quantity to diagnose the success or failure 
of explosions, because the neutrino heating rate in the 
so-called gain region, Q'^ , is proportional to (e^^) L^^ 
(e.g., equation (23) in lJankal ()2001l )). 

Figure m shows E^^^ as a function of (e^^)^Ly^. Note 
in the plot that we set the horizontal axis not as (e^^ ) Lj^^ 

but as (ei/^)^ L^^ so that we can deduce the following de- 
pendence more clearly and easiljQ- In this figure, let us 
first focus on red pluses, green crosses, and blue squares 

(•^D (=/o°°*'^x4^/>'(«i'x)//o°°*i'xeiJi„/>'(<:i'x)) and (e^^) 



can be simply connected as (e^ ) 
spectrum of equation JT} . 



2.1 (e^^) for the neutrino 
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Figure 6. Diagnostic energies for exploding models at several 
hundred seconds after the bounce. Red points, green crosses, and 
blue squares correspond to models with is =100, 150, and 200 ms, 
respectively. Red circles represent the result of 2D simulation (see 
text for details). 

whose difference is characterized by tg (2D results (filled 
circles) will be mentioned in the later section). Red 
{ts =100 ms), green {ts =150 ms), and blue {ts =200 
ms) points have a clear correlation with (ej^^ ) ^ L^^ . Or- 
ange and light-blue regions represent the non-exploding 
regions for red and blue points, respectively. Both of 
them show that the minimum E"^^^ decreases with ts, 

indicating that the critical values of {ty)^ L^^ for explo- 
sion sharply depends on ts- This is because the mass 
outside the shock wave gets smaller with time so that 
the minimum energy to blow up star gets smaller too. 
By the same reason, i?diag becomes larger as ts becomes 
smaller given the same {e^J)^ L^^. To obtain a larger 
^diat" earlier spectral swapping is more preferential. 

Figure [7] shows the neutrino heating rate and the den- 
sity distribution of NH13R30E13T100S for 10 ms and 
250 ms after tg (=100 ms after the bounce). As the 
shock wave propagates outward, the density in the gain 
region sharply drops (e.g., 100-200km, dashed blue line), 
leading to the suppression of the heating rate (dashed 
red line). This is the reason of the saturation in E^diag as 
shown in Figure ID 

The remnant mass is an important indicator to diag- 
nose the consequences of the explosion in producing ei- 
ther a neutron star or a black hole. The last two lines 
in Table [T] show the integrated masses in the regions of 
p > 10^° g cm~'^ at i = is and t = oo. The latter one is 
estimated by the fitting as 



Mio(0 = A/i?(l + e-='+'^), 



(11) 



where c and d are the fitting parameters. For the 
exploding models, becomes generally smaller 

than MIq*^ because of the mass ejection. Exceptions 
are weakly exploding models (NH13R20E15T150S, 
NH13R20E15T200S, NH13R30E13T100S, and 
NH13R50E11T100S), in which the mass accretion 
continues after ts and stops eventually at late time 
(maximum masses arc presented in Table [T]). For 
the nonexploding models, the remnant mass simply 
increases with time. Regarding the 13 AIq progenitor 
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p for 250 ms after t^ 



1 10 

lo'- 

i 10' 




2 
1.5 
1 

0.5 

X -0.5 
O 

-1 
-1.5 

-2 



Radius [km] 

Figure 7. Heating rate at 10 ms (red-solid line) and 250 ms (red- 
dashed line) after the bounce for the model NH13R30E13T100S 
and the density profile (blue-solid and dashed lines). As the density 
decreases due to the neutrino driven wind, the heating rate also 
decreases. 

ivestigated in this section, the remnant masses in models 
that produce strong explosion {E^^^ > 10^^ erg), are 
considerably smaller (1.1-1.2 Mq) if compared to the 
typical mass as of observ ed neutron stars ^ IAMq 
(jLattimer fc Prakash|[2007l ). This may simply reflect the 
light iron core (~ I.26M0) inherent to the progenitor 
or the existence of mass accretion induced by the 
matter fallback after the explosion. Now we move on 
to investigate the progenitor dependence in the next 
section. 

3.1.3. The progenitor dependence 
In addition to the 13 Mq progenitor by 



iNomoto fc Hashimoto! 



we arc going to investigate 



the progenitor dependence in ID simulations. The com- 
puted models are NH15 flS M^l (iNomoto fc Hashimotol 
mEB), sl5s7b2 (ISMq) (IWooslev fc Weaved IIQqC 
sl5.0 (ISMfT,) s2q .O (2OM0), and s25.0 (25Mq) 
(jWooslev et al.ll2002[) . which are listed in Tabled The 
first sets of characters for these models indicate the 
progenitors as, 



• NH: (jNomoto fc Hashimotol 1 1 9881) 

• WW: (IWooslev fc Weaved[T995l) 

• WHW: (jWooslev et aLllMil 

Figure [8] depicts density profiles of these progenitors 
100 ms after the bounce as a function of the enclosed 
mass {Mr). It can be seen that the density profiles for 
Mr < 0.8Mq are almost insensitive to the progenitor 
masses despite the difference in the pre-collapse phase 
(see, e.g.. Figure 1 of lBurrows et a l. 2007t|). On the other 
hand, the profiles of Mr > 0.8 A/© differ between pro- 
genitors so that the critical heating rates and E^^^ are 
expected to be different also. In Figure [HI the envelope 
of WHW25 is shown to be thickest, while the envelope 
of NH13 is thinnest. 

Figure d shows the critical heating rates as a function 
of the progenitor masses. In agreement with intuition, 
the critical heating rate for models WHW25 and NH13 
belongs to the high and low ends, respectively. However, 



N13 (kriss) 

N15 (NH88) 
S15 (WW95) 
H15 (WHW02) 
H20 (WHW02) 
H25 (WHW02) 




1 1.5 

Enclosed Mass [M,: ] 



Figure 8. Density profiles of investigated progenitors 100 ms after 
the bounce as functions of the enclosed mass. 
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Figure 9. The critical heating rate, (eiy^) Lv^, as a func- 
tion of t he pro genitor mass, M . Circles are progeni- 
tors from INomot o & Has himotol I I198SI) , the s quare is from 
IWooslev Weaveil tl995( ). and crosses are from IWooslev et al.l 
II2OO2I) . respectively. The error bars represent the distance between 
the last failing and the first exploding model in our grid of models. 
The symbols locate at the centers of error bars. The error bar is 
small for model NH13 because wc calculated a more refined grid 
of models for 13A/q progenitor (Table [TJ than for the 15-25Mq 
progenitors (Table [2]|. 

the critical heating rate for model WHW20 is almost the 
same as the one for model NH13 although the envelope 
of model WHW20 is much thicker than model NH13 (see 
Figure [8]). Our results show that the critical heating rate 
is indeed affected by the envelope mass, however, the 
relation is not one-to-one. It is also interesting to note 
that the critical heating rates for ISM© progenitors of 
WW15, WHW15 and NH15, are different by a factor of 
~ 3, which may send us a clear message that the accurate 
knowledge of supernova progenitors is also pivotal to pin 
down the supernova mechanism. 

The integrated masses with p > 10^° g cm~'^ for t = tg 
and t = 00 are listed in the last two lines in Table [2] and 
Figure [TOl The tendencies are the same as found with 
NH13. As for model WHW25, we obtain results with 
^dia<^ > 10^^ erg and M^ > 1.4Af0, simultaneously. 

3.2. Two-dimensional models 
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Table 1 

ID simulations 



Model 


Dimension 


Ru 






ts 


Explosion 


diag 

[10^1 erg] 










[km] 


[MeV] 


[lO^^gj-g s-l] 


[ms] 




[Me] 




NH13R10E15T100S 


ID 


10 


15MeV 


0.29 


100 


No 


— 


1.18 


— 


NH13R10E17T100S 


ID 


10 


17MeV 


0.48 


100 


No 


— 


1.18 


— 


NH13R10E18T100S 


ID 


10 


18MeV 


0.60 


100 


No 


— 


1.18 


— 


NH13R10E19T100S 


ID 


10 


19MeV 


0.75 


100 


Yes 


1.00 


1.18 


1.14 


NH13R10E20i 100b 


ID 


10 


20JVleV 


0.92 


100 


Yes 


1.49 


1.18 


1.12 


NH13R20E13T100S 


ID 


20 


ISJVleV 


0.66 


100 


No 


— 


1.18 


— 


NH13R20E13T150S 


ID 


20 


13MeV 


0.66 


150 


No 


— 


1.21 


— 


NH13R20E13T200S 


ID 


20 


13MeV 


0.66 


200 


No 


— 


1.25 


— 


NH13R20E14T100S 


ID 


20 


14MeV 


0.88 


100 


No 


— 


1.18 


— 


NH13R20E14T150S 


ID 


20 


14MeV 


0.88 


150 


No 


— 


1.21 


— 


NH13R20E14T200S 


ID 


20 


14MeV 


0.88 


200 


No 


— 


1.25 


— 


NH13R20E15T100S 


ID 


20 


15MeV 


1.16 


100 


Yes 


0.97 


1.18 


1.15 


NH13R20E15T150S 


ID 


20 


15MeV 


1.16 


150 


Yes 


0.54 


1.21 


< 1.24 


NH13R20E15T200S 


ID 


20 


15MeV 


1.16 


200 


Yes 


0.47 


1.25 


< 1.26 


NH13R20E21T100S 


ID 


20 


21MeV 


4.47 


100 


Yes 


5.56 


1.18 


1.07 


NH13R20E22T100S 


ID 


20 


22MeV 


5.39 


100 


Yes 


6.50 


1.18 


1.07 


NH13R28E13T100S 


ID 


28 


13MeV 


1.29 


100 


No 


— 


1.18 


— 


NH13R29E13T100S 


ID 


29 


13MeV 


1.38 


100 


No 


— 


1.18 


— 


NH13R30E11T100S 


ID 


30 


llMeV 


0.76 


100 


No 


— 


1.18 


— 


NH13R30E11T150S 


ID 


30 


llMeV 


0.76 


150 


No 


— 


1.21 


— 


NH13R30E11T200S 


ID 


30 


llMeV 


0.76 


200 


No 


— 


1.25 


— 


NH13R30E12T100S 


ID 


30 


12MeV 


1.07 


100 


No 


— 


1.18 


— 


NH13R30Eli i IoOd 


ID 


30 


12MeV 


1.07 


150 


No 




1.21 




NH13R30E12T200S 


ID 


30 


12MeV 


1.07 


200 


No 


— 


1.25 


— 


NH13R30E13T100S 


ID 


30 


13MeV 


1.48 


100 


Yes 


0.85 


1.18 


< 1.19 


NH13R30E13T150S 


ID 


30 


13MeV 


1.48 


150 


No 




1.21 




NH13R30E13T200S 


ID 


30 


13MeV 


1.48 


200 


No 




1.25 




NH13R30E14T100S 


ID 


30 


14MeV 


1.99 


100 


Yes 


1.58 


1.18 


1.12 


NH13R30E14T150S 


ID 


30 


14MeV 


1.99 


150 


Yes 


0.98 


1.21 


1.19 


NH13R30E14T200S 


ID 


30 


14MeV 


1.99 


200 


Yes 


0.68 


1.25 


1.22 


NH13R30E15T100S 


ID 


30 


15MeV 


2.62 


100 


Yes 


2.27 


1.18 


1.10 


NH13R30E15T150S 


ID 


30 


15MeV 


2.62 


150 


Yes 


1.43 


1.21 


1.16 


NH13R30E15T200S 


ID 


30 


15MeV 


2.62 


200 


Yes 


0.93 


1.25 


1.22 


NH13R30E20T100S 


ID 


30 


20MeV 


8.28 


100 


Yes 


6.84 


1.18 


1.07 


NH13R40E11T100S 


ID 


40 


llMeV 


1.35 


100 


No 




1.18 




NH13R50E11T100S 


ID 


50 


llMeV 


2.10 


100 


Yes 


0.86 


1.18 


< 1.18 


NH13R60E11T100S 


ID 


60 


llMeV 


3.03 


100 


Yes 


1.48 


1.18 


1.12 
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Figure 10. The final NS masses as a fu nction of the progeni- 
tor m ass, M. Circles a re progenitors fromlNomoto &: Hashimotol 
||1988 | ). squares are f r om | \yooslev fc Weaveil 1T995I ). and crosses are 
from fWooslev et al.l 1)20021 ). respectively. 

Here we discuss the effects of spectral swapping in 2D 
(axisymmctric) simulations. Since our 2D simulations, 
albeit utilizing the IDSA scheme, are still computation- 
ally expensive, it is not practicable to perform a system- 
atic survey in 2D as we have done in ID simulations. 
Looking at Figure [9] again, we choose models WHW15 



("Woo slev et"alll2002[ ) and NH13 (jNomoto fc Hashimoto! 
1988), whose critical heating rate belong to the high and 
low ends, respectively. 

3.2.1. 2D without spectral swapping 

The basic hydrodynamic picture is the same with ID 
before the shock-stall (e.g., till < 10 ms after bounce). 
After that, convection as well as SASI sets in between 
the stalled shock and the gain radius, which leads to 
the neutrin o -heat ed shock revival for model NH13 (e.g., 
ISuwa et all ([Mol) ). While for model WHW15, the po- 
sition of the stalled shock, following several oscillations, 
begins to shrink at > 400 ms after bounce. 

Even after the shock revival, it should be emphasized 
that the shock prop agation for mode l NH13 is the so- 
called "passive" one (jBuras "eFaII[2006l) . This means that 
the amount of the mass ejection is smaller than the accre- 
tion in the post-shock region of the expanding shock (see 
moti ons of mass s hells in the post-shock region of Figure 
1 in ISuwa et al.l (j2010l )). Some regions have a positive 
local energy (Eq. ©), but the volume integrated value 
is quite as small as < 10^*^ erg at the maximum. In order 
to reverse the passive shock into an active one it is most 
important to energize the explosion in some way. Using 
these two progenitors that produce a very weak explo- 
sion (model NH13) and do not show even a shock revival 
(model WHW15), we hope to explore how the dynamics 
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Table 2 

Progenitor dependence 



Model 


Dimension 


Rv 








Explosion 




71 — 








[km] 


[MeV] 


[lO^^erg s-i] 


[ms] 




[10^1 erg] 


[Me] 


[Mq] 


NH15R30E11T100S 


ID 


30 


llMeV 


0.76 


100 


No 




1.34 




NH15R30E12T100S 


ID 


30 


12MeV 


1.07 


100 


No 




1.34 




NH15R30E13T100S 


ID 


30 


13MeV 


1.48 


100 


Yes 


0.65 


1.34 


< 1.38 


iNxiiO-txoUrji^ 1 luuo 




Qn 
oU 


i4nvie V 


1 OQ 

i.yy 


lUU 


I es 


Z. i / 


i.o4 


1 oc; 
1 . ZO 


JNH15K30ijl5 i 100b 


ID 


30 


15MeV 


2.62 


100 


Yes 


3.73 


1.34 


1.21 


WW15R30E11T100S 


ID 


30 


llMeV 


0.76 


100 


No 




1.40 




WW15R30E12T100S 


ID 


30 


12MeV 


1.07 


100 


No 




1.40 




W W15R30E13T100S 


ID 


30 


13MeV 


1.48 


100 


No 




1.40 




VV VV LOiXoyjiLiL^ ± lUUO 




OU 


i^ivie V 


1 QQ 

I.yy 


lUU 


X es 


1 Qzl 




l.Ol 


WW15R3UE15i 100b 


ID 


30 


15MeV 


2.62 


100 


Yes 


3.41 


1.40 


1.25 


WHW15R30E11T100S 


ID 


30 


llMeV 


0.76 


100 


No 


— 


1.49 


— 


WHW15R30E12T100S 


ID 


30 


12MeV 


1.07 


100 


No 




1.49 




WHW15R30E13T100S 


ID 


30 


13MeV 


1.48 


100 


No 




1.49 




WHW15R30E14T100S 


ID 


30 


14MeV 


1.99 


100 


No 




1.49 




WHW15R30E15T100S 


ID 


30 


15MeV 


2.62 


100 


Yes 


3.55 


1.49 


1.36 


WHW20R30E11T100S 


ID 


30 


llMeV 


0.76 


100 


No 




1.45 




WHW20R30E12T100S 


ID 


30 


12MeV 


1.07 


100 


No 




1.45 




WHW20R30E13T100S 


ID 


30 


13MeV 


1.48 


100 


Yes 


0.99 


1.45 




WHW20R30E14T100S 


ID 


30 


14MeV 


1.99 


100 


Yes 


2.20 


1.45 


1.34 


WHW20R30E15T100S 


ID 


30 


15MeV 


2.62 


100 


Yes 


3.61 


1.45 


1.29 


WHW25R30E12T100S 


ID 


30 


12MeV 


1.07 


100 


No 




1.69 




WHW25R30E13T100S 


ID 


30 


13MeV 


1.48 


100 


No 




1.69 




WHW25R30E14T100S 


ID 


30 


14MeV 


1.99 


100 


No 




1.69 




WHW25R30E15T100S 


ID 


30 


15MeV 


2.62 


100 


Yes 


0.73 


1.69 


< 2.00 


WHW25R30E16T100S 


ID 


30 


16MeV 


3.39 


100 


Yes 


5.92 


1.69 


1.49 


WHW25R30E17T100S 


ID 


30 


17MeV 


4.32 


100 


Yes 


9.21 


1.69 


1.41 



would cliangc when tlic spectral swapping is switched on. 

3.2.2. 2D with spectral swapping 

Table [3] shows a summary for our 2D models, in which 
the last character of each model (A) indicates "Axisym- 
metry". Models NH13A and WHW15A are 2D models 
without spectral swapping for NH13 and WHW15, re- 
spectively. 

As in ID, the onset of the spectral swapping is taken 
to be ts = 100 ms after bounce. At this time, model 
NH13 shows the onset of the gradual shock expansion 
with a small diagnostic energy of i^diag ^ 3 x 10^^ erg, 
and the shock radius is located at ^ 300 km. As for 
model WHW15, there is no region with a positive local 
energy (e.g., Eq. ([9|)) and the shock radius is ^ 200 km. 
The density profile for this model is essentially same as 
the one in the ID counterpart (see Figure [5]) but with 
small angular density modulations due to convection. 

In Figure [6l red filled circles represent ^^^iag model 
NH13. It can be seen that the critical heating rate to 
obtain E^^^ ~ 10^^ erg is smaller for 2D than the corre- 
sponding ID counterparts (compare the heating rates for 
{e^y - 2.2 X lO^"* MeV^ erg s'^). In fact, models 
with {e^y L^^ < 2.2 x lO^"* MeV^ erg s"! fail to explode 
in ID, but succeed in 2D (albeit with a relatively small 
-^diag 1^^^ than 10^^ erg). As opposed to ID, it is rather 
difficult in 2D to determine a critical heating rate due to 
the stochastic nature of the explosion triggered by SASI 
and convection. In our limited set of 2D models, the crit- 
ical heating rate is expected to be close to {e^^ L^^^ ^ 
1.5 X lO^"' MeV^ erg s~^, below which the shock does not 
revive (e.g., {e^y L^j^ ^ 1 x 10^** MeV^ erg s~^, is the 
lowest end in the horizontal axis in the figure). 

As seen from Figure [H -Edfag becomes visibly larger 



for 2D than ID especially for a smaller (e,y^) L^_^. As 
the heating rates become larger, the difference between 
ID and 2D becomes smaller because the shock revival 
occurs almost in a spherically symmetric way (before 
SASI and convection develop non-linearly) . In Table 3, 
it is interesting to note that model NII13R3GE11T100A 
fails to explode, while we observed the shock-revival for 
the corresponding model without the spectral swapping 
(model NH13A). This is because the heating rate of 
model NH13R30E11T100A is smaller than NH13A due 
to the small (e,y^), which can make it more difficult to 
trigger Vx explosions. On the other hand, if the energy 
gain due to the swap is high enough (i.e., for models with 
greater than E12 in Table [3]), the swap can facilitate ex- 
plosions. 

Figure [TT] depicts the entropy distributions for mod- 
els NH13A (top panel) and NH13R30E13T100A (bot- 
tom panel). It can be seen that m odel NH13A show s 
a unipolar-like explosion (see also iSuwa et al.l 120101) . 
while model NII13R30E13A explodes rather in a spher- 
ical manner as mentioned above. Model NH13A expe- 
riences several oscillations aided by SASI and convec- 
tion before explosion, while the stalled shock for model 
NII13R30E13T100A, turns into expansion shortly after 
the onset of the spectral swapping. In fact, the shapes 
of hot bubbles behind the expanding shock are shown to 
be barely changing with time (bottom panel), which in- 
dicates a quasi-homologous expansion of material behind 
the revived shock. 

Figure [12] shows the time evolution of mass 
shells for models NH13A (thin-gray lines) and 
NH13R30E13T100A (thin-orange fines). Black and 
red thick lines represent the shock position at the north 
pole for each models. The mass shells for model NII13A 
continue to accrete to the PNS, since the shock passively 
expands outwards as already mentioned. Due to this 
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Figure 11. Time evolution of the entropy distributions. Top: NH13A without the spectral swapping for 100, 200, 300, and 450 ms after 
bounce from left to right. Bottom: NH13R30E13T100A for 100, 150, 200, 250 ms after bounce (corresponding to 0, 50, 100, 150 ms after 
the onset of the spectral swapping.) 

Here let us discuss a validity of the parameters for the 
spectral swap that we have assumed so far. For example, 
the criteria of explosion for model NH13R30E12T100A 
was L^^ « 1.07X lO^^erg s'^ and (e^J « 12 MeV. These 
values are even smaller than the typi cal values obtained 
in ID Boltzmann simulations (e.g., iLiebendorfer et al.l 

(|2004[ )). which show L^^ w 2x 10^^ gj-g g-i and ^^(e^J « 

20 MeV (i.e. (ci/^) ~ 14 MeV with a vanishing chemical 
potential) earlier in the postbounce phase. Therefore the 
spectral swapping, if it would work as we have assumed, 



continuing mass accretion, the remnant for this model 
would be a black hole instead of a neutron star. On 
the other hand, model NH13R30E13T100A shows a 
mass ejection with a definite outgoing momentum in 
the postshock region so that the remnant could be 
a neutron star. Unfortunately however, wc cannot 
predict the final outcome due to the limited simulation 
time. A long- t erm s imulation recently done in ID (e.g., 
iFischer et al.l (|201CI[ )) should be indispensable also for 
our 2D case. This is, however, beyond the scope of this 
paper. 



10 



SUWA ET AL. 




50 100 150 200 250 300 350 400 
Time after boune [ms] 



Figure 12. Time evolution of mass shells for NH13A (thin-gray 
lines) and NH13R30E13T100A (thin-orange lines). Blaek and red 
thiek lines represent the shock position at the north pole. 

may be a potential to assist explosions. 

It should be noted that the critical heating rate in this 
study might be too small due to the approximation of 
the light-bulb scheme. In this scheme, we can include the 
geometrical effect of the finite size of the neturinosphere 
as in Eq. ([7]), but can not include the back reaction 
by the matter, i.e. the absorption of neutrino. Some 
fraction of neutrinos, in fact, are absorbed in the gain 
region and the neutrino luminosity decreases with the 
radius. We omit this effect in this study so that the 
heating rate might be overestimated in the simulation 
with the spectral swapping. Thus, the fully consistent 
simulation including the spectral swapping is necessary 
for more realistic critical heating rate, which is beyond 
the scope of this study. 

Finally we discuss the 15 Mq progenitor labeled by 
WHW15. As mentioned, this progenitor fails to explode 
without spectral swapping even in 2E0. Figure [13] shows 
the entropy distributions of WIIW15A (left; nonexplod- 
ing) and WHW15R30E15T100A (right; exploding) for 
220 ms after the bounce (corresponding to 120 ms after 
ts for model WHW15R30E15T100A). The model with 
R^^ — 30 km and (e^^) = 14 McV does not explode 
in ID but explodes in 2D (compare Table [2] and |3]). 
Again, the mulitidimensionality helps the onset of ex- 
plosion. The critical heating rate in 2D is in the range 
of 2.5 < {e^^f L^J {10^^ MeV2 erg s'^) < 3.9, while 
it is 3.9 < {eu^ L^J(IQ'^'^ MeV^ erg s'^) < 5.9 in 
ID. Therefore the critical heating rate in 2D can be by 
a factor ^ 2 smaller than in ID. In 2D, a critical I'x 
luminosity and average energy to obtain explosion are 
^ 2x 10^^ erg s~-^ and {e^J ^ 14 MeV (correspond- 
ing to \JJ/i^ ^ 20 MeV) , which are close to the results 
obtained in a ID Boltzmann simulation (jSumivoshi et ahl 
120051 ) for a 15 Mq progenitoiQ- The diagnostic energy 

^ This is eonsist e nt w ith a very reeent result by 
lObergaulinger &: Jankal II2011I ). who performed 2D simulations of 
model WHW15 with spectral neutrino transport 

Note that the progenitor employed in lSumivoshi et al.l 120051 ) 
is WW95, so that the direct comparison may not be fair. However, 
the critical heating rate in ID for WW15 is smaller than WHW15 
(Figure[9]l and the mass of the envelope is thicker for WHW15 than 
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Figure 13. The entropy distributions of WHW15A (left) and 
WHW15R30E15T100A (right) for 220 ms after the bounce. 

as well as the estimated remnant masses are listed in the 
last three columns in Tabled E^^^ (as well as M^^^) for 
exploding models is shown to be larger than the model 
series of NII13. As a result, some of the 2D models for 
WHW15 produce strong explosions {E^^^ ~ 10^^ erg), 
while simultaneously leaving behind a remnant of 1.34- 
1.52 Mq. We think that it is only a solution acciden- 
tally found by our parametric explosion models. However 
again, the critical heating rates that require to assist the 
neutrino-driven explosion via the spectral swapping are 
never far away from the ones obtained in the Boltzmann 
simulations. We hope that our exploratory results may 
give a momentum to supernova theorists to elucidate the 
effects of collective neutrino oscillations in a more con- 
sistent manner. 

4. SUMMARY AND DISCUSSION 

We performed a series of one- and two-dimensional 
hydrodynamic simulations of core-collapse supernovae 
with spectral neutrino transport via the IDSA scheme. 
To model the spectral swapping which is one of the 
possible outcomes of the collective neutrino oscilla- 
tions, we parametrized the onset time when the spec- 
tral swap begins, the radius where the spectral swap 
takes place, and the threshold energy above which the 
spectral interchange between heavy-lepton neutrinos and 
electron/anti-electron neutrinos occurs. By doing so, we 
systematically studied the shock evolution and the mat- 
ter ejection due to the neutrino heating enhanced by 
spectral swapping. We also investigated the progenitor 
dependence using a suite of progenitor models (13, 15, 
20, and 25 Mq). With these computations, we found 
that there is a critical heating rate induced by the spec- 
tral swapping to trigger explosions, which differs between 
the progenitors. The critical heating rate is generally 
smaller for 2D than ID due to the multidimensionality 
that enhances the neu trino heating efficiency (see also 
iJanka fc Muelleilll996l ). The remnant masses can be de- 

WW15 (Figure|8ll. This indicates that our discussion above seems 
to be quite valid, although we really need ID results for WHW15 
to draw a more solid conclusion. 
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Table 3 

2D simulations 



Model 


Dimension 




(^I'x ) 




ts 


Explosion 


diag 










[km] 


[MeV] 


[lO^^erg s-i] 


[ms] 




[lO-'^ erg] 


\Mq] 


[Mq] 


NH13A 


2D 










Yes 


~ 0.1 (oscillating) 






NH13R30E11T100A 


2D 


30 


llMeV 


0.76 


100 


No 




1.18 




NH13R30E12T100A 


2D 


30 


12MeV 


1.07 


100 


Yes 


0.45 


1.18 


< 1.23 


NH13R30E13T100A 


2D 


30 


13MeV 


1.48 


100 


Yes 


1.03 


1.18 


< 1.18 


NH13R30E15T100A 


2D 


30 


15MeV 


2.62 


100 


Yes 


2.33 


1.18 


1.10 


WHW15A 


2D 










No 








WHW15R30E13T100A 


2D 


30 


13MeV 


1.48 


100 


No 








WHW15R30E14T100A 


2D 


30 


14MeV 


1.99 


100 


Yes 


1.96 


1.48 


< 1.52 


WHW15R30E15T100A 


2D 


30 


15MeV 


2.62 


100 


Yes 


3.79 


1.48 


1.34 



termined by the mass ejection driven by the neutrino 
heating, which range in I.I-I.SA/q depending on the pro- 
genitors. For our 2D model of the ISMq progenitor, we 
found a set of the parameters that produces an explosion 
with a canonical supernova energy close to 10^^ erg and 
at the same time leaves behind a remnant mass close to 
~ IAMq. Our results suggest that collective neutrino os- 
cillations have the potential to solve the supernova prob- 
lem if they occurs. These effects should be explored in 
a more self-consistent manner in hydrodynamic simula- 
tions. 

Here it should be noted that the simulations in this 
paper are only a very first step towards more realistic su- 
pernova modeling. For the neutrino transfer, we omitted 
the cooling of heavy lepton neutrinos and the inelastic 
neutrino scattering by electrons. These omissions lead to 
an overestimation of the diagnostic energy and also they 
should relax the criteria for explosion. The ray-by-ray 
approximation may lead to an overestimation of the di- 
rectional dependence of the neutrino anisotropies. A fuU- 
angle transport will give us a more correct answer (see 
lOtt et al.|[2008t iBrandt et all 1201 1[ ). IMoreover, due to 
the coordinate symmetry axis, the SASI develops prefer- 
entially along the axis; it could thus provide a more favor- 
able condition for the explosion. As several expl oratory 
simulations have been done recently (e.g., J waka mi et al.l 
[200llScheidegger et al.l[20M INordhauset ar.li20ia . 3D 
supernova models are indeed necessary also to pin down 
the outcomes of the spectral swapping. 

Finally we briefly discuss whether the oscillation pa- 
rameters taken in this paper are really valid in views 
of recent work whose focus is on clarifying the still- 
veil ed nature of collec tive neutrino oscillations. Follow- 
ing |Dua^^F^L| (|201Cl[ ). there are at least two conditions 
for the onset of collective neutrino oscillations in the case 
of inverted neutrino mass hierarchy. 

The first criteria should be satisfied in the so-called 
bipolar regime of the collective oscillation. In the regime, 
the neutrino number density should exceed the critical 
value, 



1 



(yiT^- 1)' \/2Gf (ep, 
X 



:1.4 X lO^^cm"^ 



15 McV\ 



, (12) 



where x is the fractional excess of neutrinos over antineu- 
trinos, Am^ is the characteristic mass-squared splitting 
(a typical value of ~ 2.4 x 10~'^ eV^ is employed here). 




150 200 250 
Time after bounce [ms] 



400 



Figure 14. Time evolution of x = Fi- 
number flux of i/;. 



. /Fi^s — 1, where F^. is the 



and Gf is Fermi coupling constant. By using our simula- 
tion results, we can estimate x which is often treated as 
a parameter (typically ^0.01-0.25) so far. The fol lowing 
estimation is given in lEsteban-Pretel et al.l (|2007D . that 
is X — ^i/e / — 1 in the case of vanishing F^^ , where F^. 
is the number flux of Vi. From Figure [T4l it can be seen 
that X ^ 0.2-0.3 for 100-400 ms after bounce. Since the 
typical number density in the post-shock region (r ^200- 
300 km) can be estimated as, 

Li,. 



47rr^c (ep^) 
1.1 X 10^1 cm-^ 
15 MeV\ 



1052 erg g- 



100 km 



(13) 



therefore, the first condition is satisfiecQ- 

The second criteria is related to the decoherence of col- 
lective oscillations by matter. In order to overwhelm the 
suppression by the decohenrence, the following condition 
should be satisfied 

«P,,crit ne, (14) 
where is the number density of electrons where the 

Even if x is as small as x 0.01 due to the inclusion of Vx, 
the criteria could be marginally satisfied. 
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decoherence takes place. This is equivalent to, 

ipe.crit ^ Ye- (15) 

In our ID simulation, Yp^ - (0.1 - 0.2) x for 100 km 
^ ^ ''sh, where Tgh is the shock radiutQ- Since this con- 
dition is barely satisfied, the collective oscillations in re- 
ality could modify the spectrum to some extent between 
heavy- lepton neutrinos and electron/anti-electron neu- 
trinos, however the full swapping ass umed in this study 
may be exaggerated. Verv recentl-v Fl . IChakrabortv et al.l 
()2011ai rbl) pointed out that the matter effect could fully 
suppress the spectral swapping in the accretion phase us- 
in g ID neutrino- r adiati on hydrodynamic simulation data 
of iFischer et al.l ()2010[ ). However, the current under- 
standing of the collective oscillation is not completed 
and calculations in this field employ several assump- 
tions (e.g., single an gle approximation) (but see also 
iDasgupta et al.ll2011l for more recent work). To draw 
a robust conclusion, one needs a more detailed study 



including the collective neutrino flavor oscillation to the 
hydrodynamic simulations in a more self-consistent man- 
ner, which we are going to challenge as a sequel of this 
study. 
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APPENDIX 
CODE VALIDITY 
Conservation of Energy 

In this section, we demonstrate the conservation of physical quantities using the spherical collapse model (NH13). 
Figure flSl depicts the evolution of total binding energy by gravity (red line), total internal energy (green), total kinetic 
energy (blue), total trapped- neutrino energy (magenta), total energy leaked by neutrinos (cyan), and variation of 
overall energy (black dashed) , respectively. The gravitational energy and total energy are negative and absolute values 
are shown. The gravitational energy and internal energy dominate (with different sign) and reach 10'^^ erg soon 
after bounce. Despite such an enormous energy change, the total energy varies only within '--^ 3 x 10^^ erg so that the 
violation of energy conservation remains < 0.03%. The energy of the trapped neutrinos decreases with the diffusion 
timescale, which leads to the PNS cooling. The kinetic energy rapidly drops because of the photodissociation of iron 
and the electron capture {I'e emission) that is consistent with the shock stall. We have monitored these values in a 2D 
simulation and obtained a similar level of energy conservation. 




0.1 0.2 0.3 0.4 0.5 

Time after Bounse [s] 



Gravitational Binding Energy Trapped-Neutrino Energy 

Internal Energy Leaked Energy by Neutrinos 

Kinetic Energy Variation of Total Energy 

Figure 15. Time evolution of gravitational energy (red), internal energy (green line), kinetic energy (blue line), trapped-noutrino energy 
(magenta line), released energy by neutrinos (cyan line), and summation of these energies (black dashed line). These quantities are 
determined by integration with respect to volume included in our simulation except for released energy by neutrinos (magenta), which is 
f{Ljj^ + Li,^)dt. Since gravitational energy and total energy are negative, the absolute values are shown. The violation of total energy 
(dashed line) remains < 3 X 10*^ erg, which is ~ 0.03% of gravitational energy and internal energy after bounce (~ lO^"^ erg). 



Outside the shock, Yy^ > Ye is achieved due to rapid density 
decrease. 

In fact they posted their papers on astro-ph after our submis- 
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Comparison with AGILE 

Here, we present the result of our n umerical simulation in spherical symmetry and compare with the result of AGILE- 
IDS A code (jLiebendorfer et ani2009[ ). AGILE (Adaptive Grid with Implicit Leap Extrapolation) is an implicit general 
rclativistic hydrodynamics code that evolves the Einstein equations based on conservative finite differencing on an 
adaptive grid. We employ a one-dimensional version of our ZEUS-2D code that has been developed to perform 
multidimensional supernova simulations. 

We compare the evolution of a 13 Af© star of iNomoto fc Hashimotol (jl988l ) in Newtonian gravity from precollapse 
model to 100 ms after bounce. We find good agreement between the results of the ZEUS-2D and AGILE during the 
early postbounce phase when the neutrino burst is launched and the accretion shock expands to its maximum radius. 
The hydrodynamic quantities are shown in following figures. 




Figure 16. Density (left top), entropy (right top), velocity (left bottom), and electron fraction (right bottom) as a function of enclosed 
mass for the result of ZEUS-2D (red lines) and AGILE (green lines). The comparison is shown at the time just after the bounce. A 
difference is seen in the entropy profile, which comes from the difference of shock capturing scheme. 
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